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摘 要 : 饱和 水 汽 压 差 是 土壤 -植被 -大 气 连续 体 水 分 传输 过 程 的 关键 影响 因素 ,在 全 球 气候 变化 背景 下 , 预测 西 
北 地 区 饱和 水 汽 压 差 ,对 于 植被 恢复 和 农林 业 气 象 灾害 风险 评估 具有 重要 的 现实 意义 。 基 于 西北 五 省 (区 )1990 一 


2019 年 月 饱和 水 汽 压 差 值 ,采用 趋势 分 析 和 人 小波 分 析 等 方法 研究 了 西北 地 区 饱和 水 汽 压 差 年 际 变化 特征 和 周期 性 


变化 规律 ;采用 指数 模型 和 ARIMA 模型 ,筛选 最 佳 样本 步 长 和 预测 步 长 ,对 西北 地 区 饱和 水 汽 压 差 进行 模拟 和 预 
测 。 结 果 表 明 :(1) 西北 五 省 (区 ) 中 ,新疆 年 均 饱 和 水 汽 压 差 最 高 ,其 次 为 宁夏 .陕西 .甘肃 和 青海 ; 近 30 a 整 体 上 西 
北 地 区 饱和 水 汽 压 差 呈 上升 趋 势 , 其 中 宁夏 和 新 疆 饱 和 水 汽 压 差 上 升幅 度 最 大 ,分 别 为 0.036 kPa. (10a)! FI 0.033 


kPa.(10a) ,其 次 为 甘肃 [0.026 kPa-(10a)'] .青海 [0.021 


kPa- (10a)'] 和 陕西 [0.012 kPa(10a)']; (2) 西北 各 省 (区 )， 


16 a 尺度 周期 对 小 波 方 差 贡 献 最 大 ,为 饱和 水 汽 压 差 变 化 的 主 周期 。 此 外 ,陕西 .甘肃 和 新 疆 还 存在 24~27 a 的 周期 
特征 ,方差 贡献 较 小 ;(3) 相对 于 指数 模型 ,ARIMA 模型 均 方 根 误差 平均 减少 42.3% ,决定 系数 尼 平 均 提高 11.1%， 
Nash-Sutclife 效 率 系数 平均 提高 17.7% ,有 效 提 高 了 饱和 水 汽 压 差 预测 精度 ;(4) 未 来 一 段 时 间 内 ,西北 各 地 区 饱和 水 
汽 压 差 均 存在 不 同 程度 的 升 高 趋势 ,以 宁夏 和 新 疆 地 区 的 饱和 水 汽 压 差 增幅 最 为 明显 ,分别 为 9.5% 和 8.9%。 
关键 词 : 饱和 水 汽 压 差 ， 指 数 平滑 ARIMA 模型 ， 预 测 ; 西北 地 区 


饱和 水 汽 压 差 是 量化 大 气 十 燥 程 度 的 重要 指 
标 ,是 陆地 生态 系统 植物 群落 发 挥 生 态 功能 的 主要 
了 驱动 力 ,也 是 干旱 引起 植物 死亡 的 重要 驱动 因素 "。 
随 着 气温 升 高 ,饱和 水 汽 压 显著 增 大 ,水 汽 压 差 也 
增 大 2 。 饱 和 水 汽 压 差 通过 改变 植物 的 气孔 行为 从 

影响 蒸腾 ,在 一 定 条 件 下 饱和 水 汽 压 差 驱动 植被 
蒸腾 作用 ,超过 一 定 阔 值 ,会 导致 植物 气孔 导 度 下 
降 和 光合 作用 降低 ,甚至 造成 农作物 减产 和 林木 
死亡 “1。 人 研究 表明 : 近 几 十 年 来 ,全 球 饱 和 水 汽 压 
差 急 剧 增长 ,这 将 进一步 改变 土壤 -大 气 系统 水 
分 分 配 ”, 加 剧 土壤 温度 升 高 和 土壤 水 分 散失 ,增加 
了 植物 干旱 胁迫 的 风险 "1。 

我 国 西北 地 区 位 于 亚 欧 大 陆 内 地 ,远离 海洋 ， 
降水 稀少 ,植被 覆盖 率 低 , 属 典型 的 气候 敏感 带 和 
生态 脆弱 区 吓 。 自 20 世 纪 90 年 代 起 ,西北 地 区 升 
温 明 显 22 ,潜在 蔡 散 量 增加 品 。 近 几 十 年 来 ,西北 
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全 区 降水 量 呈 微弱 的 上 升 趋势 ,而 东部 区 域 降水 量 
明显 下 降 , 干 时 发生 频率 增加 ”。 全 球 气候 变化 
下 , 单 以 温度 和 降水 难以 清晰 地 描述 区 域 自 然 气候 
变异 ,而 饱和 水 汽 压 差 能 够 综合 反映 大 气 的 温度 和 
湿度 状况 ,指示 地 区 气候 的 干燥 程度 。 同 时 ,饱和 
水 汽 压 差 也 是 影响 净 生 态 系统 交换 量 ” ERA 
量 “ 以 及 植物 叶片 气孔 变化 ”的 主要 因素 。 因 此 ， 
饱和 水 汽 压 差 的 变化 ,不 仪 会 影响 陆地 植被 生长 ， 
而 且 可 能 会 对 区 域 水 循环 产生 重要 影响 。 众 多 人 研 
究 关 注 于 气候 变化 条 件 下 的 植被 动态 "生态 系统 
碳 交 换 “ 和 蒙 散发 ”特征 ,以 及 植被 蒸腾 和 气孔 导 
度 对 饱和 水 汽 压 差 变化 的 响应 ,对 于 我 国 西北 
地 区 的 饱和 水 汽 压 差 特征 .模拟 和 预测 研究 鲜 有 报 
道 。 因 此 ,选择 适用 性 好 的 模型 对 其 进行 模拟 和 预 
测 ,无 论 是 对 于 未 来 气候 变化 情境 下 植被 格局 演变 
研究 ,还 是 农林 业 气 象 灾害 风险 评估 都 尤为 重要 。 
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传统 的 经 验 回归 模型 需要 输入 的 参数 较 多 , 且 
所 需 数据 量 很 大 ,计算 过 程 繁 杂 ;神经 网 络 在 气 
候 预 测 方面 具有 一 定 的 优势 ”, 但 是 输出 结果 对 参 
数 的 敏感 性 高 ,而 部 分 参数 多 依靠 经 验 公式 确定 ， 
这 极 大 的 影响 了 模型 的 预测 能 力 ;灰色 系统 进行 气 
候 预 测 不 失 为 一 种 好 方法 ” ,但 对 于 振荡 序列 误差 
较 大 ” ,在 一 定 程度 上 限制 了 其 气候 模拟 和 预测 方 
面 的 应 用 。 指 数 模型 是 预测 时 序 未 来 值 的 常用 模 
型 ,结构 简单 ,短期 预测 能 力 好 ; 自 回归 求 积 移动 平 
均 模 型 (Autoregressive Integrated Moving Average 
model,ARIMA) 能 根据 非 平 稳 时 序 历史 数据 对 未 来 
值 预测 , 且 不 依赖 外 部 变量 ,能 够 在 一 定 程度 上 克 
服 外 部 参数 引起 模型 精度 不 足 的 问题 ” ,在 粮食 产 
量 预测 7 干旱 预测 下 和 温度 预测 ”等 方面 均 取 得 
了 较 好 的 预测 结果 。 本 研究 以 西北 五 省 (区 ) 
1990 一 2019 年 月 饱和 水 汽 压 差 数据 为 样本 ,分 析 西 
北 地 区 饱和 水 汽 压 差 年 际 变化 特征 和 周期 性 变化 
规律 ;引入 指数 模型 和 ARIMA 模型 ,通过 比较 不 同 
模型 结构 ,建立 适宜 模拟 西北 各 地 区 饱和 水 汽 压 差 
的 最 优 模型 ;通过 比较 不 同样 本 长 度 和 预测 长 度 对 
模型 精度 的 影响 ,筛选 出 预测 模型 的 最 佳 样本 步 长 
和 预测 步 长 ,提高 模型 的 精度 ,并 对 未 来 一 段 时 间 
饱和 水 汽 压 差 进行 预测 ,以 期 为 未 来 应 对 气候 变化 
条 件 下 的 干旱 半 干 旱地 区 生态 建设 和 农林 产业 合 
理 布局 提供 科学 依据 。 


1 研究 区 概况 与 研究 方法 


1.1 研究 区 概况 

西北 地 区 包括 陕西 (31°42'~39°35'N, 105°29'~ 
111°15'E) HIN (32°11'~42°57'N ,92°13'~108°46’E) , 
宁夏 回族 自治 区 (35°14’ ~39° 14’ N, 104° 17’ ~ 
109°39'E) | FY (31°36’ ~39° 19’ N, 89°35'~103°04'E) 
Ail Br E AE OK AA X (34° 25’ ~48° 10’ N, 73°40’ ~ 
96°18'E), FRE Bs AMERRE SE BE 
究 区 属于 干旱 半 干 旱地 区 ,以 温带 大 陆 性 气候 为 
主 ,年 均 温 8 % ,年 均 降 水 量 400 mm 以 下 ,由 东 到 西 
递减 。 主 要 地 形 以 山地 ` 盆 地 和 高 原 为 主 , 植 被 覆 
HAD , Fe LAY A) Ae BRE AS a 
12 数据 与 模型 构建 
12.1 数据 来 源 通过 国家 气象 数据 中 心 获取 
1990 一 2019 年 各 省 (区 ) 共 计 118 个 气象 站 点 的 温度 


和 相对 湿度 月 观测 数据 。 各 省 实际 月 饱和 水 汽 压 
差 (VPD ) 计 算 公 式 如 下 5 


VPD=0.61 te osl - RH) (1) 

式 中 :7 为 大 气温 度 (%);RH 为 相对 湿度 (% ) 。 
1.2.2 验证 序列 平稳 性 和 确定 d 值 ”饱和 水 汽 压 差 
是 温度 和 相对 湿度 的 函数 , 受 季节 等 因素 影响 呈 周 
期 性 变化 ,具有 一 定 的 趋势 项 和 季节 效应 。 因 此 初 
步 采用 Holt 指 数 平 滑 吕 和 Holt-Winters 指数 平滑 六 
对 样本 数据 进行 模拟 ;建立 ARIMA(p,d,g) 模 型 前 ， 
需 判 断 时 序 的 平稳 性 ,主要 判断 方法 有 自 相 关 图 观 
察 和 单位 根 检 验 。 其 中 自 相 关 图 观察 法 虽然 操作 
简单 ,但 不 够 严格 ,因此 ,本 文采 用 单位 根 检验 法 对 
序列 稳定 性 进行 检验 9, 以 确保 验证 的 可 靠 性 ,并 
判断 序列 数据 达到 平稳 所 需 的 差分 值 。 对 于 非 稳 
定 序列 数据 需要 进行 差分 处 理 , 即 确定 d 值 。 
1.2.3 确定 p 值 和 g 值 p 和 gq 分 别 为 ARIMA(p,d， 
dg) 模型 中 的 自 回归 项 数 和 移动 平均 项 数 。 为 筋 选 最 
TRE BAG ,本 文采 用 判断 4 次 差分 后 平稳 序列 自 
相关 (ACF) 和 偏 自 相关 (PACF) 系 数 的 截 尾 、 拖 尾 和 
周期 性 ,并 结合 建立 不 同 参 数 模型 进行 比较 的 方法 
来 确定 最 佳 模 型 结构 。 
1.3 模型 评价 指标 

调整 的 尼 表 征 了 模型 的 拟 合 优 度 ,因此 在 模型 
选择 时 采用 调整 的 尼 作 为 参考 依据 。 同 时 结合 最 
小 信息 量 AIC (Akaike Information Criterion ) 准则 来 
确定 最 优 模 型 ,调整 的 尼 越 大 ,AIC 越 小 ,模型 效果 
越 好 ;采用 统计 量 Ljung-Box Q 来 检验 残 差 序 列 是 否 
显著 自 相 关 , 若 已 > 0.05 ,说 明 拟 合 后 的 残 差 序 列 为 
和 白 品 声 序列 ;采用 绝对 误差 绝对 值 平均 (AE) 、 相 对 
误差 绝对 值 平均 (RE) 、 均 方 根 误差 (RMSE ) 评 价 模 
型 参数 ;采用 Nash-Sutclife 效 率 系数 (NSE ) 对 模型 进 
行 评 价 。 其 中 NSE 被 用 来 验证 实测 值 与 模拟 值 间 
的 拟 合 情 况 ,该 值 越 接近 于 1 ,模型 的 效率 越 高 , 适 
用 性 越 强 。 各 指标 公式 如 下 : 


AIC =-21nL+ 2k (2) 
AE=|ET, -ET „| (3) 
RE=|(ET, - ET,)/ET, x 100% (4) 
ET. -ET,,) 
RMSE = DET -ET (5) 
N, 
ET.- ET) 
NSE = 1.0- LET.- ETy) (6) 


2 


X ET, E ET ive) 
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SUP : LARRY A WSR KR; k Fy BE FY th SB 
BX N JIPE; ET, A BE ; ET wy KWE ; ET ave 
为 ET 的 平均 值 。 
1.4 数据 处 理 

采用 Excel 2010 软件 对 数据 进行 预 处 理 ; 采 用 
SPSS 19.0 和 Eviews 7.2 软件 进行 数据 分 析 和 建 模 ; 
采用 Origin 9.4 软件 绘图 ;采用 LSD (Least Significant 
Difference) 法 对 饱和 水 汽 压 差 值 进行 多 重 比 较 ; 采 
用 小 波 方差 确定 序列 变化 的 主 周 期 。 


2 结果 与 分 析 


2.1 西北 地 区 饱和 水 汽 压 差 变化 特征 

由 图 1 可 以 看 出 ,1990 一 2019 年 ,各 省 (区 ) 年 均 
饱和 水 汽 压 差 为 新 疆 (XJ) > 宁夏 (NX) > 陕西 
(SN) > 甘肃 (GS) > 青海 (QH)。QH 的 饱和 水 汽 压 
差 年 际 波动 相对 较 小 ,而 SN GS, NX 和 XJ 四 省 (区 ) 
的 饱和 水 汽 压 差 年 际 波动 较 大 。 各 省 (区 ) 饱 和 水 
汽 压 差 年 际 变 率 分 别 为 0.32%(SN) .0.24%(GS) 、 
0.32% (NX) 、0.19% (QH) 和 0.29% (XJ) 。 对 各 省 
(区 ) 年 际 饱和 水 汽 压 差 进行 线性 趋势 判断 ,NX 和 
XJ 年 际 倾向 斜率 较 高 ,分别 为 0.036 kPa (10a) ' (hs) 
和 0.033 kPa*(10a)"'(k;)。 其 次 为 GS、QH 和 SN ,分 别 
为 0.026 kPa (10a) ' (k2) 、0.021 kPa (10a)! (ha) 和 
0.012 kPa- (10a) "(ki)。 说 明 30 a 来 宁夏 和 新 疆 的 饮 
和 水 汽 压 差 整 体 水 平 上 升 趋势 比 其 他 省 份 明显 ,其 
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注 :~hs 分 别 表示 陕西 .甘肃 .宁夏 .青海 和 新 疆 的 他 和 水 汽 压 差 倾 
问 斜 率 ; 字 母 SN 表示 陕西 省 ,GS 表示 甘肃 省 ,NX 表示 宁夏 回族 自 
治 区 ,QH 表示 青海 省 ,XJ 表示 新 疆 维吾尔 族 自治 区 。 下 同 。 
1 1990 一 2019 年 西北 五 省 (区 ) 饱 和 水 汽 压 差 
Fig. 1 VPD of five provinces (regions) in Northwest 
China from 1990-2019 


他 各 省 的 饱和 水 汽 压 差 上 升 趋势 较 小 。 

各 省 (区 ) 人 饱和 水 汽 压 差 演变 过 程 存在 周期 性 
振荡 特征 (图 2) ,表现 为 年 代 际 的 变化 , 且 各 省 (区 ) 
的 周期 性 变化 规律 较为 相似 ,主要 表现 为 15~20 a 
的 振荡 周期 ,此 周期 的 振荡 周期 信号 能 力 最 强 , 周 
期 性 变化 最 为 明显 。 该 尺度 周期 在 整个 研究 期 表 
现 稳 定 且 具有 全 域 性 ,说 明 近 30 a 来 西北 地 区 饱和 
水 汽 压 差 周期 性 变化 表现 较为 平稳 。 小 波 方差 可 
以 反映 饱和 水 汽 压 差 序 列 的 波动 随时 间 尺 度 的 分 
布 ,可 用 来 确定 饱和 水 汽 压 差 变 化 过 程 存在 的 主 周 
期 。 由 图 2 还 可 以 看 出 ,各 省 (区 ) 在 5a、16 a 时 间 尺 
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图 2 各 省 (区 ) 饱 和 水 汽 压 差 小 波 分 析 


Fig. 2 Wavelet coefficients and wavelet variance of VPD among the five provinces (regions) 
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度 存 在 极 值 ,16 a 左右 的 尺度 周期 对 小 波 方差 贡献 £ 0.8 
率 最 高 , 达 7.75% ,为 饱和 水 汽 压 差 变化 的 主 周期 。 ® ，。 
同时 ,陕西 .甘肃 和 新 疆 还 存在 24-27 a 的 周期 特 £ 
征 , 但 方差 贡献 率 较 小 ,为 2.80%~3.14%。 民 04 
由 图 3 可 以 看 出 ,西北 五 省 中 ,XJ 月 最 高 饱和 水 时 
汽 压 差 达 2.34 kPa, 其 次 为 NX、SN、GS 和 QH, 分 别 “ae ee ke aa Se 
y 1.68 1.42 .1.15 kPa Fil 0.78 kPa。 各 省 (区 ) 年 均 饱 省 份 


和 水 汽 压 差 以 XJ 最 高 (图 4) ,为 0.61 kPa, 其 次 为 
NX, SN, GS 和 QH, 分 别 为 0.54、0.48、0.46 kPa 和 
0.36 kPa。 由 图 5 可 以 看 出 ,1 月 和 2 月 SN 月 均 饱 和 
水 汽 压 差 无 显著 差异 ,其 余 各 省 (区 ) 各 月 均 存在 显 
著 性 差异 (P< 0.05) ,7 一 12 月 各 省 (区 ) 饱 和 水 汽 压 
差 各 月 均 存 在 显著 性 差异 (P<0.05)。12、1 月 和 2 
月 GS 和 NX 无 显著 性 差异 ,11、12、1 月 和 2 月 SN 和 
CS 存在 显著 性 差异 (P<0.05)。 另 外 ,1 月 和 2 月 GS 
和 NX 无 显著 性 差异 ,11、12、1 月 和 2 月 SN 和 GS 存 
在 显著 性 差异 (P< 0.05)。 
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省 份 
图 3 1990 一 2019 年 西北 五 省 (区 ) 月 饱和 水 汽 压 差 特 征 
Fig. 3 Characteristics of monthly VPD in five provinces 
(regions) of Northwest China from 1990-2019 


图 4 1990 一 2019 年 西北 五 省 (区 ) 年 饱和 水 汽 压 差 特征 
Fig.4 Characteristics of annual VPD in five provinces 
(regions) of Northwest China from 1990-2019 


2.2 模型 选择 及 参数 率 定 

由 表 1 可 以 看 出 ,QH 和 XJ 饱和 水 汽 压 差 原 始 
序列 单位 根 统 计 (7) 均 大 于 1%、5% 和 10% 水 平 的 
值 ,因此 不 能 拒绝 原 假设 , 即 序 列 非 平稳 。NX 和 GS 
饱和 水 汽 压 差 原始 序列 了 值 小 于 5% 和 10% 水 平 的 
值 ,大 于 1% 水 平 的 值 ,说 明 在 1% 水 平 下 不 能 拒绝 原 
假设 。SN 饱 和 水 汽 压 差 原 始 序 列 7 值 均 小 于 1%、 
5% 和 10% 水 平 的 值 , 但 仍 存在 0.1% 的 概率 接受 原 
假设 。 对 序列 一 阶 差分 后 ,各 省 (区 ) 饱 和 水 汽 压 差 
原始 序列 7 值 均 小 于 1%、5% 和 10% 水 平 的 值 ,接受 
原 假设 概率 为 0, 说 明 在 1%、5% 和 10% 水 平 下 拒绝 
原 假设 。 因 此 可 以 确定 序列 平稳 ,ARIMA 模型 参数 
d 取 值 为 1。 

由 图 6 可 以 看 出 ,一 阶 差分 后 各 省 平稳 序列 自 
相关 系数 (ACF) 均 在 1 期 后 衰减 为 小 值 波动 的 过 程 
较为 突然 ,可 视 为 一 阶 截 尾 ,由 此 可 以 确定 一 阶 差 
分 序列 服从 MA(1)。 对 于 SN、GS、NX 和 QH, 偏 自 
相关 系数 (PACF) 在 1 期 后 衰减 波动 较 大 ,3 期 后 逐 


25 
2.0 SN ZA GS SY NX BY QH E45 XJ Af 
Š i Ae Bf 四 
Bd Ad Cd e 
C Df 
e 1.0 ABb Be Ac Ac Ce Cd D E D ! Ef 
sp 0.5} AaBaBaBaca A2BbBCbCbpp ZN 
Z 00 LE as FT 
ie 20 1 Aa 2 Ab 3 4 5 6 
g is CaCa Ba Cha, Bb = Ac 
m™ 1.0 a i Cb Db CeCe Be CdcaBd Ad ABe BCe Af 
0.0 POR ， 
7 8 9 10 11 12 


注 : 小 写字 母 表示 同一 省 份 各 月 份 水 汽 压 差 差 异 显著 性 ,大 写字 母 表示 同一 月 份 各 省 份 


水 汽 压 差 差 异 显著 性 ,不 同 字 母 表示 差异 显著 (P<0.05)。 
图 5 1990 一 2019 年 西北 五 省 月 均 饱 和 水 汽 压 差 比较 
Fig. 5 Comparison of monthly mean VPD in five provinces of Northwest China from 1990-2019 
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Tab.1 ADF test results of VPD series among the five provinces (regions) 


307 


Pa 


SN GS NX QH XJ 
原始 序列 T -4.20 -2.87 -3.01 -2.52 -2.05 
1% -3.45 -3.45 -3.45 -3.45 -3.45 
5% -2.86 -2.86 -2.86 -2.86 -2.86 
10% -2.57 -2.57 -2.57 -2.57 -2.57 
AP/% 0.10 4.90 3.50 11.20 26.56 
一 阶 差 分 T -17.68 -15.18 -15.36 -18.36 -22.37 
1% -2.58 -2.58 -2.58 -2.58 -2.58 
5% -1.94 -1.94 -1.94 -1.94 -1.94 
10% -1.62 -1.62 -1.62 -1.62 -1.62 
AP/% 0.00 0.00 0.00 0.00 0.00 


主 :AP 表 示 接 受 概率 ;SN 表示 陕西 省 ,GS 表 示 甘 肃 省 ,NX 表 示 宁 夏 


回族 自治 区 ,QH 表示 青海 省 ,XJ 表 示 新 疆 维 吾 尔 族 自治 区 。 下 同 。 


渐 在 置信 区 间 内 小 范围 波动 ,因此 可 认为 服从 AR 
(1) AR(2) 或 AR(3)。 对 于 XJ,PACF 在 2 期 后 快速 
衰减 至 置信 区 间 范 围 波动 ,因此 可 以 考虑 AR(2) 或 
AR(3)。 综 上 所 述 ,初步 选择 ARIMA 模型 如 表 2 所 
IR ,对 不 同 模型 进行 比较 ,最 终 科 选 出 SN GS FI XJ 
的 最 优 ARIMA 模型 为 :ARIMA (3,1,1) ;NX 为 ARI- 
MA(1,1,1);QH 为 ARIMA(2.1,1)。 筛 选 的 最 优 指数 
模型 如 表 3 所 示 。 
2.3 模型 模拟 及 效果 分 析 

模型 残 差 检 验 发 现 , 当 指数 模型 样本 步 长 至 少 
为 4a,ARIMA 模型 样本 步 长 至 少 为 6a 时 ,模型 才 
有 效 。 因 此 对 指数 模型 设置 了 步 长 为 4~30 a 间隔 
为 2a 的 样本 数据 ;对 ARIMA 模型 设置 了 步 长 为 6~ 
30 a 间隔 为 2a 的 样本 数据 。 经 分 析 比 较 ,筛选 出 
各 模型 最 佳 样本 步 长 ( 表 4)。 为 实现 序列 的 连续 性 
并 对 各 省 (区 )2020 年 及 以 后 的 饱和 水 汽 压 差 进行 
预测 ,因此 确定 了 指数 模型 拟 合 时 段 为 2016 一 2019 
年 ,ARIMA(3,1,1) 和 ARIMA(2,1,1) 模 型 拟 合 时 段 为 
2004—2019 年 , ARIMA (1,1,1) 模 型 拟 合 时 段 为 
1998 一 2019 年 。 

设置 预测 步 长 分 别 为 3.5. 7a 和 9a, 以 最 佳 样 
本 步 长 建立 模型 并 进行 回归 分 析 ( 表 5)。 预 测 步 长 
为 3a 时 Holt 指 数 平滑 的 尼 和 NSE 均 表现 最 高 , 随 
着 预测 步 长 增加 拟 合 度 下 降 。 不 同 预测 步 长 下 ， 
Holt-Winters 指数 平滑 模型 的 尼 相 差 不 大 ,但 在 预测 
步 长 为 3a 时 误差 最 小 , 且 模 型 的 NSE 最 大 ,说 明 模 
型 效率 更 高 。ARIMA(3,11) 和 ARIMA(1,1,1) 模 型 
分 别 在 预测 步 长 为 7a 和 5 a 时 拟 合 度 和 模型 效率 更 
高 。ARIMA(2,1,1) 模 型 在 各 预测 步 长 条 件 下 尼 无 


明显 差异 ,但 步 长 为 7a 时 误差 表现 更 小 ,NSE 更 高 。 

综 上 所 述 ,指数 模型 最 佳 预测 步 长 为 3 a; ARI- 
MA(3,1,1) 和 ARIMA(2.1,1) 模 型 最 佳 预测 步 长 为 7 
aiARIMA(1,1,1) 模 型 最 佳 预测 步 长 为 Sa。 相 对 于 
指数 模型 ,ARIMA 模型 的 RMSE 平 均 减 少 42.3% , R? 
平均 提高 11.1% ,NSE 平 均 提 高 17.7%。 
2.4 模型 预测 

根据 筛选 的 最 佳 预测 步 长 ,确定 指数 模型 的 预 
M R JE X 2020—2022 4E ; ARIMA (3,1,1) FI ARIMA 
(2,1,1) 模 型 的 预测 尺度 为 2020 一 2026 年 ;ARIMA 
(1,1,1) 模 型 的 预测 尺度 为 2020 一 2024 年 。 对 各 省 
(区 ) 饱 和 水 汽 压 差 进行 预测 ,结果 显示 ( 表 6) ,未 来 
3 a(2020 一 2022 年 ) 陕 西 和 新 疆 的 饱和 水 汽 压 差 将 
分 别 增 加 3.6% 和 6.7%, 而 甘肃 .宁夏 和 青海 的 饱和 
水 汽 压 差 将 有 所 减 小 ;未 来 5a(2020 一 2024 年 ), 宁 
夏 饱 和 水 汽 压 差 增加 9.5%;2020 一 2026 年 ,陕西 . 甘 
肃 .青海 和 新 疆 饱 和 水 汽 压 差 均 呈 不 同 程度 的 升 高 
趋势 ,分 别 增 加 3.6% 、8.0%、7.9% 和 8.9%。 


3 讨论 


饱和 水 汽 压 差 是 用 于 模拟 植被 气孔 时 度 ” 和 
佑 算 蒸 散 量 ? 等 生态 模型 的 重要 气象 参数 之 一 ,也 
是 生态 系统 碳 通 量 及 水 汽 通 量 研 究 的 主要 参考 变 
HE 饱和 水 汽 压 差 升 高 可 能 会 限制 陆地 生态 系 
统 中 植物 碳 吸收 和 水 分 利用 。 部 分 研究 表明 ,近年 
来 全 球 饱 和 水 汽 压 差 急 剧 增长 ,可 能 是 由 于 海洋 燕 
发 减少 ,陆地 水 汽 含量 降低 ”。 全 球 气候 变 暖 背景 
下 ,近年 来 西北 地 区 气温 升 高 ,极端 高 温 事件 发 生 
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图 6 饱和 水 汽 压 差 一 阶 差 分 序列 自 相 关 
Fig. 6 First-order differential sequence autocorrelation 
diagram of VPD 


频次 增加 ,这 将 可 能 引起 饱和 水 汽 压 差 增 大 。 
本 研究 结果 显示 ,1990 一 2019 年 ,宁夏 和 新 疆 的 饮 
和 水 汽 压 差 升 高 趋势 明显 ,这 可 能 与 地 区 温度 升 
高 ,干燥 事件 发 生 频次 增加 有 关 ”"。 降 水 主要 通 
过 影响 空气 温度 和 湿度 来 影响 饱和 水 汽 压 差 中 ,而 
饱和 水 汽 压 差 的 变化 会 影响 陆 面 实际 蒸 散发 史 , 从 

影响 区 域 的 降水 状况 。 近 年 来 ,我 国 西北 东部 区 
域 降 水 有 明显 下 降 趋势 中 ,说 明 降水 的 减少 在 一 定 


表 2 ARIMA 模型 效果 对 比 
Tab.2 Comparison of ARIMA model effects 


省 份 模型 JH R AIC 
SN ARIMA(1,1,1) 0.45 -0.32 
ARIMA(2,1,1) 0.46 0.25 
ARIMA(3,1,1) 0.47 0.43 
GS ARIMA(1,1,1) 0.39 -0.78 
ARIMA(2,1,1) 0.45 -0.02 
ARIMA(3,1,1) 0.46 0.22 
NX ARIMA(1,1,1) 0.38 -0.04 
ARIMA(2,1,1) 0.38 0.70 
ARIMA(3,1,1) 0.37 0.92 
QH ARIMA(1,1,1) 0.43 -1.80 
ARIMA(2,1,1) 0.47 -0.99 
ARIMA(3,1,1) 0.47 -0.76 
XJ ARIMA(2,1,1) 0.49 1.81 
ARIMA(3,1,1) 0.50 2.06 


表 3 指数 模型 效果 对 比 


Tab.3 Comparison of exponential model effects 
省 份 指数 模型 调整 的 尼 均 方 根 误差 
SN Holt 0.67 0.130 
GS Holt-Winters 0.70 0.081 
NX Holt-Winters 0.69 0.120 
QH Holt-Winters 0.71 0.046 
XJ Holt-Winters 0.68 0.103 


表 4 模型 最 佳 样本 步 长 
Tab.4 Optimum sample step among the five models 


模型 最 佳 样 本 步 长 /a 均 方 根 误差 平均 绝对 误差 


Holt 4 0.086 0.067 
Holt-Winters 4 0.076 0.053 
ARIMA(3,1,1) 16 0.113 0.079 
ARIMA(2,1,1) 16 0.048 0.036 
ARIMA(1,1,1) 22 0.130 0.092 


程度 上 对 饱和 水 汽 压 差 变 化 起 到 了 正 向 作用 ,使 得 
饱和 水 汽 压 差 增长 趋势 明显 。 本 文 对 西北 各 省 
(区 ) 饱 和 水 汽 压 差 进行 预测 ,研究 表明 ,2020 一 
2026 年 各 省 (区 ) 饱 和 水 汽 压 差 均 呈 上 升 趋势 ,这 可 
能 意味 着 西北 地 区 的 暖 干 化 趋势 仍 会 维持 较 长 一 
段 时 间 。 有 研究 表明 ,在 温度 与 降水 的 双重 作用 
下 ,整个 西北 地 区 仍 呈 变 干 的 趋势 ”, 这 与 本 研究 
结果 基本 一 致 。 另 外 ,陕西 饱和 水 汽 压 差 年 际 上 升 
趋势 较 小 ,而 年 际 变 率 相对 较 大 ,因此 短期 的 预测 
可 能 具有 不 确定 性 。 但 对 各 省 (区 ) 的 预测 结果 进 
行 综 合 评 价 分 析 , 从 饱和 水 汽 压 差 与 温度 变化 直接 
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表 5 不 同 预 测 步 长 下 模型 输出 参数 比较 
Tab.S Comparison of output parameters of models with different predicted steps 

模型 模拟 步 长 /a R 均 方 根 误差 相对 误差 /% 绝对 误差 /kPa NSE 
Holt 3 0.76 0.1809 21.3 0.1196 0.6713 
5 0.75 0.1935 22.3 0.1327 0.6383 
7 0.73 0.1924 22.5 0.1320 0.6233 
9 0.73 0.1910 23.2 0.1347 0.6292 
Holt-Winters 3 0.95 0.0939 9.8 0.0584 0.9207 
3 0.96 0.1003 11.7 0.0688 0.9128 
7 0.95 0.1080 12.9 0.0734 0.9013 
9 0.94 0.1154 14.4 0.0802 0.8901 
ARIMA (3,1,1) 3 0.97 0.0926 17.1 0.0735 0.9543 
5 0.96 0.0878 15.9 0.0699 0.9566 
7 0.97 0.0854 17.7 0.0697 0.9589 
9 0.96 0.0878 17.4 0.0718 0.9553 
ARIMA(2,1,1) 3 0.94 0.0537 9.6 0.0363 0.9096 
5 0.94 0.0512 9.2 0.0363 0.9114 
7 0.94 0.0461 8.7 0.0326 0.9256 
9 0.92 0.0497 8.8 0.0347 0.9154 
ARIMA(1,1,1) 3 0.92 0.1111 15.0 0.0852 0.9098 
3 0.93 0.1063 14.9 0.0794 0.9250 
7 0.93 0.1081 14.9 0.0807 0.9243 
9 0.92 0.1024 14.7 0.0806 0.9222 

表 6 指数 模型 和 ARIMA 模型 的 预测 结果 
Tab. 6 Prediction results of exponential models and ARIMA models 
SN 尺度 /a GS 尺度 /a NX 尺度 /a QH 尺度 /a XJ 尺度 /a 

MV/(kPa: mon”) 0.55 = 0.50 - 0.63 - 0.38 - 0.90 - 

PV-E/(kPa: mon") 0.57 3 0.43 3 0.61 3 0.32 3 0.96 3 

PV-A /kPa.mon ) 0.57 7 0.54 7 0.69 5 0.41 7 0.98 T 


注 :MV 表示 实测 值 ;PV-E 表 示 指 数 模型 预测 值 ;PV-A 表 示 ARIMA 模 型 预测 值 。 下 同 。 


相关 这 一 角度 出 发 ,这 与 前 人 的 研究 结果 具有 较 好 
的 一 致 性 中 。 以 2016 一 2019 年 实测 饱和 水 汽 压 差 
和 预测 饱和 水 汽 压 差 平均 水 平 进 行 比较 ( 表 7) ,发 
现 指数 模型 预测 值 在 整体 上 相对 于 ARIMA 模型 偏 省 份 NS 0 NW N 
低 , 即 指数 模型 预测 结果 会 低估 饱和 水 汽 压 差 的 平 tes > ie = 
均 水 平 。 ARIMA 模型 预测 结 采 更 接近 于 实际 水 平 ， PV-A /(kPa: mon) 0.56 0.54 0.68 0.41 0.95 
却 并 不 能 反映 出 饱和 水 汽 压 差 短 时 间 内 的 波动 变 一 
化 ,但 分 别 从 指数 模型 和 ARIMA 模 型 对 新 疆 地 区 饮 
和 水 汽 压 差 的 预测 结果 来 看 , 均 高 于 过 去 阶段 的 饱 
和 水 汽 压 差 水 平 ,预计 在 未 来 一 段 时 间 内 (2020 一 


表 7 水 汽 压 差 预 测 值 与 实测 值 水 平 比较 


Tab.7 Comparison between predicted values and 


measured values of VPD 


资源 合理 分 配 等 来 降低 这 种 风险 。 
指数 平滑 和 ARIMA 模型 广泛 应 用 于 各 个 人 研究 
领域 ”1 ,尤其 是 ARIMA 模 型 的 预测 能 力 具 有 显 车 


2026 年 ) ,新 疆 地 区 的 饱和 水 汽 压 差 可 能 持续 上 
升 。 而 长 期 保持 这 种 水 汽 压 差 升 高 趋势 ,将 可 能 

生 潜 在 的 土壤 干旱 和 植被 退化 的 风险 ,因此 蝇 需 采 
取 相 应 的 措施 ,如 通过 农林 产业 结构 调整 ,进行 水 


优势 ””。 李 俐 等 ”采用 ARIMA 模型 和 径 向 基 神 经 
网 络 预测 夏 玉米 单产 量 , 结 果 表 明 ARIMA 模型 精度 
Eo Duan 等 等 采用 指数 平滑 模型 预测 地 下 水 
位 取得 较 好 结果 。 本 文 应 用 指数 模型 和 ARIMA 模 
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型 对 西北 五 省 (区 ) 饱 和 水 汽 压 差 进行 模拟 和 预测 ， 
结果 表明 ,相对 于 指数 模型 ,ARIMA 模型 均 方 根 误 
差 平均 减少 42.3% ,决定 系数 尼 平 均 提 高 11.19 ， 
Nash-Sutclife 效 率 系数 平均 提高 17.7%,ARIMA 模型 
的 拟 合 度 更 高 ,预测 能 力 更 好 。 

影响 饱和 水 汽 压 差 变化 的 因素 错综复杂 ,包括 
与 其 直接 相关 的 大 气温 度 和 湿度 ,而 大 气温 度 和 湿 
度 又 受 经 纬度 海拔、 太阳 辐射 以 及 人 类 活动 等 多 
种 要 素 的 影响 。 饱 和 水 汽 压 差 预测 精度 的 关键 是 
模型 选择 ,为 进一步 提高 模型 的 适用 性 和 预测 
度 , 在 今后 的 研究 中 可 尝试 更 多 的 优化 模型 和 组 合 
模型 ,以 期 增强 模型 对 饱和 水 汽 压 差 趋 热 变 化 的 识 
别 能 力 。 另 外 ,样本 和 预测 步 长 均 会 对 预测 结果 产 
生 影 响 。 本 文通 过 对 比分 析 不 同 模型 .不 同样 本 步 
长 和 预测 步 长 对 模型 产生 的 误差 ,筛选 最 优 模型 结 
构 .最 佳 样本 步 长 和 预测 步 长 ,建立 的 饱和 水 汽 压 
差 指数 预测 模型 可 解释 76% 以 上 的 变量 ,纳什 效率 
系数 达 0.6 以 上 ;ARIMA 模型 可 解释 93% 以 上 的 变 
量 , 纳 什 效 率 系数 达 0.9 以 上 ,有 效 提 高 了 饱和 水 汽 
压 差 的 预测 精度 。ARIMA 模型 结构 简洁 直观 ,准确 
率 高 ,适用 性 好 ,可 为 西北 地 区 饱和 水 汽 压 差 的 预 
测 提供 参考 。 


y 


ah 


4 结论 


(1) 近 30a 来 ,整体 上 西北 地 区 饱和 水 汽 压 差 
旺 上 升 趋势 。 宁 夏 和 新 疆 的 饱和 水 汽 压 差 增幅 最 
为 明显 ,分 别 为 0.036 kPa:(10a)' 和 0.033 kPa+(10a)", 
其 次 为 甘肃 .青海 和 陕西 ,分 别 为 0.026 .0.021 kPa 
(10a) FH 0.012 kPa.(10a) ;西北 各 省 (区 ) 饱 和 水 汽 
压 差 存在 以 5a 和 16 a 为 中 心 尺度 的 周期 ,16 a 周期 
对 方差 贡献 最 大 , 达 7.75%, 为 饱和 水 汽 压 差 变化 的 
主 周期 。 

(2) 根据 现 有 观测 资料 ,预计 未 来 7a(2020 一 
2026 年 ) 西 北 地 区 饱和 水 汽 压 差 继续 升 高 。 至 2026 
年 ,陕西 .甘肃 .青海 和 新 疆 的 饱和 水 汽 压 差 将 分 别 
增加 3.6% .8.0% .7.9% 和 8.9%; 至 2024 年 ,宁夏 饱和 
水 汽 压 差 增 加 9.5% ;甘肃 .宁夏 和 青海 地 区 的 饱和 
水 汽 压 差 在 2020 一 2022 年 有 所 回落 ,但 总 体 仍 呈 上 
升 趋势 。 

(3) 相 较 于 指数 模型 ARIMA 模型 误差 小 , 拟 
合 度 好 ,对 变量 的 解释 度 高 ,适用 性 强 , 能 更 好 的 对 


饱和 水 汽 压 差 进行 预测 ,在 今后 的 研究 中 ,ARIMA 
模型 可 以 作为 西北 地 区 气候 预测 模型 选择 的 重要 
参考 。 
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Prediction of vapor pressure deficit in Northwest China based on 
exponential and ARIMA models 


HAN Yonggui'’’, HAN Lei'”’, HUANG Xiaoyu GAO Yang"? 
(1. School of Recourses and Environment, Ningxia University, Yinchuan 750021, Ningxia, China; 
2. China-Arab Joint International Research Laboratory for Featured Resources and Environmental 
Governance in Arid Regions, Yinchuan 750021, Ningxia, China; 3. Institute of Environmental 


Engineering, Ningxia University, Yinchuan 750021, Ningxia, China) 


Abstract: Vapor pressure deficit (VPD) is a key factor affecting water transport in the soil- plant- atmosphere 
continuum; in the context of global climate change, predicting VPD has practical significance for vegetation 
management and risk assessment of meteorological disasters affecting agriculture and forestry in Northwest 
China. Using VPD data from five provinces (regions) in Northwest China from 1990 to 2019, we analyzed the 
characteristics of VPD interannual variation and periodic variation using trend and wavelet analyses. The optimal 
sample step and prediction step were selected; exponential models and autoregressive integrated moving average 
(ARIMA) models were used to simulate and predict VPD in Northwest China. Among the five provinces, 
Ningxia had the highest trend slope of VPD [0.036 kPa- (10a) '], followed by Xinjiang [0.033 kPa:(10a) ']. The 
annual average VPD in Xinjiang was the highsest at 0.61 kPa, followed by Ningxia, Shaanxi, Gansu, and Qinghai 
(0.54 kPa, 0.48 kPa, 0.46 kPa, and 0.36 kPa, respectively). Over the past 30 years, the VPD in Northwest China 
followed an upward trend; Ningxia and Xinjiang had the largest increases in VPD at 0.036 and 0.033 kPa:(10a) `, 
respectively, followed by Gansu [0.026 kPa-(10a) '], Qinghai [0.021 kPa:(10a) '], and Shaanxi [0.012 kPa- (10a)`]. 
Compared with the exponential model, the root mean square error (RMSE) of the ARIMA model was reduced by 
42.3%, the R’ was increased by 11.1%, and the Nash-Sutcliffe efficiency coefficient increased by 17.7%. Thus, 
the VPD prediction accuracy was effectively improved. The VPD in Northwest China is expected to increase by 
varying degrees; Ningxia and Xinjiang showed the highest VPD growth rates of 9.5% and 8.9%, respectively. 
Keywords: vapor pressure deficit; exponential smoothing; ARIMA model; prediction; Northwest China 


